o 



Riesz exponential families on 
homogeneous cones 



o 

^ ; I. Boutouria*, A. Hassairi*^ 

tí 

K-i , 

Abstract. In this paper, we introduce, for a multiplier X) a notion of generalized power 
function X ^ A^(x), defined on the homogeneous cone "P of a Vinberg algebra Á. We 
then extend to A the famous Gindikin result, that is we determine the set of multiphers 
' X such that the map 6 i— > A-^(0~^), defined on V* , is the Laplace transform of a positive 

measure R^. We also determine the set of x such that R-^ generates an exponential family, 
i-tí ■ and we calcúlate the variance function of this family. 
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^ 1 Introduction 

O , 

Q\ \ It is well known (see Casalis and Letac (1996)) that the Wisliart distributions on 

the cone of (r, r) positive symmetric matrices or on the symmetric cone f2 of any 
Euclidean Jordán algebra E of rank r are the elements of the natural exponential 
^ . families generated by the measures /ip such that the Laplace transform is defined on 

üby 

for p in {-, 1) 2' ' ' '' — 2 — ^'^"^ — 2 — ' "^^^ measure /ip is absolutely continuous 

r — 1 

when p g] , +oo[ and is singular concentrated on the boundary of the cone, 

1 3 r — 1 

when p G {-,!,-,■ ■ ■, — ^ — }. In 2001, Hassairi and Lajmi have introduced the 

Riesz distribution on f2 as an extensión of the Wishart distribution. These authors 
have started from the fact that in a Jordán algebra, besides the real power of the 
determinant, there is the so called generalized power As(x) of an element x of Í2, 
defined for a fixed ordered Jordán frame of E and for s = (si, ■ ■ ■, s^.) in IR*", and they 
have used a remarkable result, due to Gindikin (1964), which determines the set S 
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of s in ]R'' such that As{6~'^) is the Laplace transform of some positive measure Rg 
on E. The gencralizcd powcr As(x) is a power function of the principal minors of 
X which reduces to (dct(x'))'' in the particular case where Si = S2 = ■ ■ ■ — Sr = p, 
and in this case, the measure Rg in nothing but ¡ip. We mention here that Ishi 
(2000) has given a more detailcd dcscription of the Gindikin set S based on the 
orbit structure of íl under the action of somc Lie group. He has also given explicitly 
the measure R^ for cach s in S. In all thcsc works. the definition of the Ricsz measure 
Rg and in particular of the Riesz probability distribution is based on the cholee of 
a totally ordered Jordán frame which allows the definition of the principal minors 
and of the generalized power of an element of the algebra. The fact that the order 
is total is a fundamental condition not only for the definition of the distribution 
but also in the proof of many results. To define models in which some specified 
conditional independencies, usually given by a graph, are taken into account, there 
has been an interest in probability distributions on the homogeneous cone of a 
Vinberg algebra. For instance, Andersson and Wojnar (2004) have defined a class 
of absolutely continuous "Wishart" distributions on an homogeneous cone. These 
distributions have been characterized by Boutouria (2005, 2007) in the Bobecka and 
Wesolowski (2002) way. They have also been charcterized by Boutouria and Hassairi 
(2008) in the way given in Olkin and Rubin (1962) for the ordinary Wishart. The 
aim of the present work is to use an approach similar to the one used in the definition 
of a Riesz exponential family on a symmetric cone to introduce a Riesz exponcntial 
family on an homogeneous cone. The distributions in these families are defined 
for any graph, that is for any order relation not necessary total. Some of these 
distributions are absolutely continuous with respect to the Lebesgue measure and 
some are singular concentrated on the boundary of the cone. In this connection, 
the Riesz distribution on the symmetric cone of a Jordán algebra may be seen as 
the particular one corrcsponding to the particular directed graph with vértex set 
{1, • ■ ■,r} and edges defined by the usual order on integers. We first define for an 
element of an homogeneous cone two kinds of principal minors, minors which are 
said strict and minors which are said large. We then define for a multiplier x, a 
notion of generalized powcr function x 1— A^{x). One of our main results is the 
determination of the set of multipliers x such that the map 6 1— > A^(^~^) is the 
Laplace transform of a positive measure R^. It is a generalization of Gindikin result 
with a more elabórate proof adapted to the properties of the Vinberg algebra and 
the graph. Concerning the generated exponential families, we give a necessary and 
sufficient condition on x in order that R^ generates an exponential family and, under 
this condition, we determine the variance function of the family. 



2 Vinberg algebras and homogeneous cones 

In this section, we introduce some notations and review some basis concepts con- 
cerning Vinberg algebras and their homogeneous cones. We also introduce a useful 
decomposition of an element of the cone. 
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Throughout the paper, I denote a partially ordered finitc sct cquipped with 
a relation denoted ^. We will write i -< j ií i ^ j and i ^ j. For all pairs 
E I X I with j -< i, let Eij be a finite-dimensional vector space over IR with 
riij — dim{Eij) > 0. Set 

{R, for i = j 
Eij for j 'í^ i OI j -<i 
{0} otherwise. 

and A= W Aij- An element A = (a¿j, i,j e /) of ^ may be seen as a matrix 

and so we define the trace trA = '^au. We also define 

iei 

"^i- = = '^rif^i, ni = l + ^(n^. +n,i),i G / and n, = ^n,- (2.1) 

Let fij : Eij — > E^j, i y j, he involutional hnear mappings, i.e., f¡^^ = fij. They 
induce an involutional mapping { A t-^ A* ) oí A given as foUows: A* = {a*j\{i,j) e 
I X I), where 

Íau for i = j 

fijiüij) = a*j for j -<i OI i -<j 
{0} otherwise. 

Let % = {A = {üij) eA, el ^ = 0},% = {A = {üij) e A, Wi,j e 

I ■ j 7^ i ^ CLij — 0} and H — {A E A, A* — A} denote respectively the set of upper 
triangular matrices, the lower triangular matrices and the Hcrmitian matrices. The 
scts of upper and lower triangular matrices in V with positive diagonal elements are 
respectively denoted by T^^ and 7¡^. The sets of diagonal matrices and of diagonal 
matrices with positive enfries are denoted by V and P"*", respectively. 

The space A is equipped with a bilinear map called multiplication and denoted 
by {A, B) I— > AB, using bilinear mappings Aij x Ajk — >■ Aik, denoted by (ajj, hjk) i— > 
ttijbjk, such that AB — C = {cij\{i,j) e / x /) with Cjj — aiijhnj- 

The multiplication is required to satisfy the foUowing properties: 

i) 'iAeA; A ^ O ^ ii{AA*) > O 

ii) y A, B eA; (AB)* = B*A* 

iii) VA, B eA] ii{AB) = tiiBA) 

iv) yA,B,CeA; tT{A{BC)) = tT{{AB)C) 

v) yU,S,TeTi; {ST)U = S{TU) 

vi) yU,S,Terf, T{UU*) = {TU)U*. 

An algebra A with the above structurc and properties is called a Vinberg algebra 
(For more details, we can refer to Andersson and Wojnar (2004). Define the inner 
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products on Eij, i >~ j by || aij ||¿^= aijfij{aij), aij G Eij. Thus instead of 

specifying the bilinear form (a¿j, bji) i-^ üijbji on Eij one can specify an inner product 
on ii^jj, i >- j. It can be established that the foUowing two conditions also 
must hold: 

1. V ttij e Eij, bjk e Ejk : II a¿j-6jfe ||¿fc=|| aij ||¿_^-|| 6^^ lljfe, k <j -<i, 
and 

2. If Oifc e Eik, bjk e í^jfc, with k ^ j <i and (a^fc, Cijbjk)ik = O for all G Eij, 
then [duaik, cijbjk)ik — O for alH G / with i -< í, and all cij G E'íj, and dn G E'h- 

We consider the elemcnt {aij\{i,j) G / x /) of "D such that a¿j = 1, V z G / as the 
unit element of A and we denote it by e. We also define Ek = {dij)ij^i G V with 
dkk = 1 and = O Vj 7^ /c. It is clear that y^^Ek = e. 

kei 

Vinberg (1965) proved that the subset V = {TT* e A, T e Ti+} C H G A 
forms a homogeneous cone, that is the action of its automorphism group is transitive. 
Let G be the connected component of the identity in Aut(P); the group of linear 
transformations leaving V invariant. We recall that x : G 1— > IR+ is said to be a 
multipher on the group G if it is continuous, = 1 and x{9i92) = x{9i)xÍ92) for 
all gi,g2 G G. Consider the map i: : T & Tf^ ^ t^ÍT) G 7r(7^''') C G such that for 
X = VV* G 1/ G 7;+ 

n(T)(X) = (TV){V*T*). (2.2) 

Andersson and Wojnar (2004) have shown that the restriction of a multiplier x to the 
(lower) triangular group T^"*", i.e., x o tt : T^"*" — > IR+ is in one to one correspondence 
with the set of (A¿,¿ G /) G K^. We will then describe a multipher x by its 
corresponding point in and we denote X = {x '■ — ^ 

If ^ °PP is the opposite ordcring on the Índex set /, i.e., i ^ j ^ j di i- The 
Vinberg algebra ^ °pp = Aij'^ , where 

IR for i = j 

Eij for j >- °PP i OI j -< °PP i 
{0} otherwise, 

differs from the Vinberg algebra A only in the ordering of the Índex set /. Vinberg 
(1965) proved that P^opp = {T*T E A, T e 7;+} is the dual cone of P. The inner 
product {A, B) — > ti{AB) on H identifies H with its dual H* , i.e., 

H ^ H* 

A ^ (B^tr(AB)), 
and this isomorphism identifies °pp with the dual cone V* oíV. 



A°PP - ) 
■^13 — \ 
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Now, for í G /, we denote = {j El; j ^ i} and Jj^ = {j G /; i ^ j} and we 
say that j separates ii and Í2 if j G /¿i^ fl /¿j^ and j ^ {ii, Í2}- In this case, j is 
called a separator. We denote Si = {j G /¿^; j is a separator} and S = \^ Si. 

iei 

lí T = {tij)ijQi is in Ti, we define the element Tj^ of 7/ by 

Ti^ — {tij)ijQi, with t'ji^ — tjk ií i di j, k and í^-j. = O otherwise, (2.3) 
and the element Tj_. of % by 

Tj_. = it'ij)ijeii with í^j;. = tjk if i -< j, A; and t'j,^ = O otherwise. (2.4) 

líX ^ TT*, we denote 

-'^¿^ — Ti^T*^, Xi^ — Ti^T*^- (2.5) 

We also denote by Vi^ (resp Vi^) the set of = T¿^7;% (resp = T¿^7;*^) 
corresponding to T G Ti'^. It is easy to see that Vi^ and P,^ are respectively the 
homogeneous cones of the Vinberg subalgebras of A defined by Ai^ — Am and 

Ai^ = denote by e¿ and é¿ respectively, the unit element of Ai^ and 

Ai^. Wc also define the rank of Vi^ (resp the rank of Vi^) the cardinal of the set 
{j ^ I-ii d j} (resp the cardinal of the set {j E I,i ^ j}). Finally, if we denote 
p = {i E 1; li^ and I^i = 0} and if we set, for X E V, 



Xi^i 



Xi^ - J2 if « e p 

X,._ ií íes (2-6) 

o otherwise. 



then, we have the foUowing decomposition of X 

X = Y.X,. (2.7) 



3 Riesz measures on an homogeneous cone 

The definition of a Riesz measure on the symmetric cone of a Jordán algebra relies 
on the notion of generalized power of an element of the cone which is a power 
function of the so-called principal minors. In order to define a Riesz distribution on 
an homogeneous cone, we need to extend all these things to a Vinberg algebra. 



3.1 Generalized power 

Wc first introduce a notion of determinant. For X — TT*, with T — (tij) E Ti, we 
define the determinants 

detX^Yltl, det-X^i^ [] 4' det^X^i= [] íj- 
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For X — TT*, with T E %, we define the strict principal minor of order k oí X a,s 

¡ det^iX^,) if 7^,^0 

1 if j^^^0 (3.8) 

and the large principal minor of order A; of X as 

A^,(X) = dct^(X^,). (3.9) 

Definition 3.1 Let x — {\, i & 1} be a multiplier andX e V, then the map defined 
by 

X^A,(X) = n(^^)^^ (3.10) 

is called the generalized power function corresponding of x- 
We also denote 

Note tliat, if Aj = A, i G /, then A^(X) = (dctX)'^. It is easy to verify that 
^y,+AX) = Ax(^)A^,(X), where ^ + x' = {A^ + A^, i G /}• 

Example 3.1 Let us consider / = {1, 2, 3, 4} and the poset defined by 

1^3, 1^4, 2^3. 

For X = TT* e V, with T = {tij) G we have A^i(X) = det-{X^i) = 
tl^, A^2{X) = det^{X^2) = tl^, ^<a{X) = det^iX^^) = tl^tj^ and A^^iX) = 
det-{X^3) = ífiíisíia- Henee, /or x = {A¿ = 1, i G I}, 

-f2 j-2 j-2 j-2 j-2 j-2 j-2 



A ("r'i — "11 22 '^11'^ 44 '^ll'^22'^33 _ j.2 .2 .2 .2 
^Xl^^ ~ 1 1 ^2 +2 +2 — ''ll''22'^33''44. 
± 1 til ''11 ''22 



3.2 Orbit decomposition of the closure V oí V 

For i G pUíS, we denote by the set of maps í/j defined from I into {0, 1} as foUows: 

If i G p. ú' is such that = O, when j ^ /¿^ or j G 5, and if i G 5, '?/^ is such 
that = O, whcn j ^ Jj^. Similarly, we denote by the set of maps ip defined 
from / into {0, 1}. lí i & p, ip is such that = O, when j ^ /¿^ or j G 5, and 
if i G 5, 'i/^ is such that ■ip{j) = O, when j ^ U^. With these notations, we define 
for ¿ G p U <S and G £^ = á\Q.g{^p) = áia.g{^p{j), j G /) G P and we denote by 
Tj^.e^ = {Te^T*, T G T^"*"}. We also consider the two elements of A 

^i^í ei in Ai^ ¿¿ ^ í in Ai^ 

1 O elsewhere , 1 O elsewhere . 

Next, we state and prove a fundamental result. It is a decomposition of V in orbits. 
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Theorem 3.1 i) 

P= E ^í^- (3-12) 

ii) Let i & pU S, then 

"Pi^ = U (3-13) 

Proof i) Let Z eV, thcn thcrc cxist a scqucnce {Z^'^^^n&j^^ in V such that Z^^^ 
Z as n ^ oo. Since {Z^^^^neN is in using the decomposition (2.7) we write 

= J2Z¡"-\ where z}"^ e Henee Z - E^¿' ^^^^re Z¿ e and (3.12) is 

i&i iei 

proved. 

ii) We will prove (3.13) by induction on the rank of thc cone Vi^. It is obvious 
that (3.13) holds for i e p U 5 such that rankPj^ = 1. Suppose that (3.13) holds 
for any i e p U 5 such that rankPj^ < / and let us show that it holds for i such 
that rankT'j^ = l. Consider the set — {j e I; I^j — {i}}- Then using the 
decomposition defined by (2.6) and (2.7) for an element of the cone Vi^, we easily 
see 

"Pi^ = E Pj<- (3-14) 

jeMi 

As rankPj^ = /, we have that rankPj^ = 1 — 1 and it foUows that rankT^^^ < ¿ — 1, 
Vj e Mj. Using the induction hypothesis, we can write 

Pj^ = U j e Mr 

Now, let = E ^® obtain 

= E = E ( U V-e^) = U ^^-e^ C U ^^-e^- (3-15) 

To conclude, we will verify that for Z e Vi^, there exist í/j E and T e 7^"*" such 
that Z = T.e^. Let Z e T'j^, then there exists a sequence {Z'^'^^}neN in Vi^ such 
that ^ Z as n ^ oo. As e there exists [7^5^ = {uf^)j^kei in ^ such 
that = C/ÍJ(C/¿5V (s^® (2-3) and (2.5)). In particular, we have 

4Í = (4Í r + E II 4? III,, (3.16) 

for k G /j^. This implies that the sequences (■u[,'^,'')„gN and (■u^"'')„gN are bounded. 
Thereforc there exists a subsequence of positivo integers (n^) such that («¿1'"'*)™ 
and (uu'r^)m converge. Let ílhi, = lim nl'l'"-' and íti.,- = hm ■ul""''. Then 

hm f/í-"""^ = Üi^, so that Z = Üi-^Ü*^ = {zkj)kiei- As ^íjj > O, we will con- 
sider separately the case zu = O and the case zu > 0. 
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Suppose that zu = 0. Then uu = {zuY^"^ = O, so that Zki — uaUki = 0, i ^ k. Thus 
Z = Zi^ G Vi^, and the result follows according to (3.15). 
If Zu > O, we consider the elements of Ai^ 

"'=(ó< ej (i £/!)■ 

where C¿ = ^üji and ¿/¿^ = {üjkjjMh^- Then üi^ = ü^üi and {¿¿j = (^^íí)^/^ > 0. 

Let Ti and T2 in 7^^, such that Ti = -u* in Ai^ and T2 = «i in Ai^, we have 
Tj.é' e Pj^. By induction hypothesis, there exists a unique ■^i e such that 
■01 (i) = O and there exists Tj e 7^"*" such that T2.É^ — í^.e^^. Let € such that 
V'2(i) = 1 and V2(j) = O Vj 7¿ ¿ and put = V'i + "02 e and T = T1T2 e %+. 
Then we have 

Z = TiT2.{Ei + É^) 
= Ti.Ei + T2.É' 
— Ti.E'j + Tj.e^j 
= (TiT2).e^ 
= T.e^, 

and (3.13) is proved. □ 
3.3 Gamma functions 

We use the generahzed power function to introduce a generahzed gamma function 
on an homogeneous cone. 

For i e p U 5, e and Xi = {^j^ 3 ^ -^i = O, Vj ^ we set 

A'(0) = {x¿ e A" I Aj- = O, for all j e such that V(j) = 0}. (3.17) 

For every Xi £ '^(0); we define a generahzed power function on Tj'^.e^ by 

Aj^(T.e^) = Ag(TT*), yTeT+, (3.18) 

where Ag(TT*) is defined by (3.11). We also define n'^ = (nj , j e 7) by 



4 = E^(^)^fcj y^eií^. (3.19) 

When 7¿ O, we introduce the measure i^,^ on Ti'^.e^ defined by 

u^diT.e^)) ^ A'-ÜíTT*) n dtkj, (3.20) 

0(j) = 1 
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where xt = {^j ^ ^ I-, such that A,- = —i¡){j){\ + n^j)/'2, if j G /¿^ and \j = 
O if j ^ li^}, and T = {tjk)j,kei ^ Finally, we denote by uq be the Dirac measure 
at 0. 

Theorem 3.2 Leti e pUS and Xi = {^j, j e A^- O, Vj ^ e X{ip). The 
integral 

r7¡+.e,(x.)= L exp{-M^)}AÍ,(^)^^(d^) (3.21) 
converges if and only if Xi £ ^ii^) satisfies the following condition: 

Xj > ^ yj & I such that = 1. (3.22) 

Moreover, under this condition, one has 



r^-.e,(x.) = 2-l%-l"'^l/^ n r(A,-^), (3.23) 



where = ^ ipU) l^'^l = X/ ^j. 



Proof If t/) = O, the integral (3.21) reduces to 1. Thus (3.22) and (3.23) hold 
trivially. lí íp ^ O, then writing Z = U.e^, where U — {ujk)j^kei ^ the integral 
(3.21) can be written 

r-7¡+.e^(x¿) = / + exp{-( Y. («ii + II II «fej llfej))} n ^S' ^dujjdukf 
' V'(i)=i 3<k i<3 

m = 1 

For j e li^, let 

Cj = {sin)i,n&i ^C^ = Y1 ^kj, with s¿j = uij ií j -< l and Si„ = O otherwise. (3.24) 
It is clear that || Cj |p= ^ || H^^- and dCj = du^j. Henee 

V;(j) = 1 V'(j) = 1 



Therefore the convergence condition is reduced to the one corresponding to the 

2 



ordinary gamma functions, that is \j > \/j E I such that ■0(j) = 1. □ 
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Remark 3.1 From Theorem 3.2, we have, for i E I , a relation between r-j-+ ^ (xi) 
and Tplx). In fact, if we denote by li the element of e"^, such that ip{j) = 1, V j G 
/j^ \S, if i E p and such that = 1, V j e /j^, if i E S, then it is olear that 

and using (2.1), for x = Yxi ^ ^, where Xi = {Aj, j G /; A^- = O, Vj ^ li^}, we 
have 

n r^-.e.,(xO = 2-i^i7r^ nr(A. - f ) = 2-i^irp(x)- 
3.4 Riesz measures 

For the definition of the Riesz distribution, we need to introduce some other nota- 
tions. Let i e pUiS and Xi = {Aj, j G /; Xj — O, Vj ^ /j^} and introduce for ^0 e 
and a; e the foUowing sets 

V') = I Xi e A*; A, = O for j ^ h^, Xj = ^ when j e /¿^ and = o| , (3.25) 



B{i, a;) = <^ Xj e A"; v^- = O for j ^ k^^'dj = when j e /¿^ and a;(j) = O L (3.26) 



S(z, i¡J) = lx^ e i3(^, ^), A, > ^ when j e and ^(j) = 1 1 , (3.27) 



S(^, í^) = I X¿ ^(^ ^) Aj > Y i ^ ^(i) = 1 [ , (3.28) 



5(z) = U E{i, V^), S(i) = U E{1 u), 5 = E Eiz). (3.29) 

For every Xi = {Xj, j G /} G A, let 

Xi = {Xj - (1 - V'(j))^,when j e 7,^, and O if j ^ 7,^}. (3.30) 

It is clear that if Xi £ B{h i^), then Xj e A'(-í/'). 

In what foUows, we denote the Laplace transform of a positive measure ¡i on the 
cone V by 

LJe) = í exp{-tr(eZ)}ii(dZ), 6 e V*. (3.31) 

7-p 
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Theorem 3.3 There exists a positive measure such that the Laplace transform 
is defined on V* and is equal to A^{9~^) if and only if x^"^- 

The proof of Theorem 3.3 rehes on the foUowing proposition. 

Proposition 3.4 Let i e p\J S. Then there exists a positive measure R^. such 
that the Laplace transform is defined on V* and is equal to AW(^-i) if and only if 

Proof Let Xí G "(O- Then there exists í/j E such that Xi £ ip). It is clear 
that Xi defined by (3.30) satisfies (3.22). We will show that the Laplace transform 
of the measure 

RxÁdZ) = \ At(Z)l ,^^^{Z)u^{dZ) 
is defined on V* and is given by 



In fact, as 9 E V*, then defined by (2.5) is in the dual cene V*^ of Vi^, and 
there exists T in 7^^ such that 9i^ = T*.E^. Let Y = 7r(T)(Z), where tt is defined 
by (2.2). As Z e 7;+.e^, there exists 5 G 7;+ such that Z = S.e^. This with (3.20) 
imply that 

p^{dY) = p^{d{n(T){Z))) = A^l{9)p^{dZ), (3.32) 

where xf = {Aj G IR, j G / such that = (1 — '4>{j))-Yi if i ^ and O if j ^ -^i^}- 
Since Xi £ ^{i^)-, then 



^ í 7\ — 

= Aj.(r-^ev,)Aj.(l^) 



A|(Z) = A|(7r-i(r)(F)) 



= A^^(r^)AÍ (y). (3.33) 
Using (3.30), (3.32) and (3.33), we get 

At{Y)v^[dY) = A^^{9)At{Z)v^[dZ). (3.34) 

Then 

Lr^^O) = exp{-tr(^7r-^(r)(y))}Ag(r^)A^.(n^V>(^n 

= Ag(r^)^-J^ /^^^^ exp(-try}Aí (y).,(.y) 
= Ag(r^). 
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Suppose that there exists a positivo measure R^^ such that the Laplace transform 
is dcfincd on V* and is equal to A^*](í^^^). Our aim to show that Xi ^ '^i'^)- 
For this Xi s-íid a t/' in consider the generahzed positive Riesz measure which 
we also denote R^^ defined for ¡f in the Schwartz space S{A) oí rapidly decreasing 
functions on A by 

RxÁf)- j. \^, L ^{Z)At{Z)u4dZ). (3.35) 

Wc will prove by induction on the rank of the cone Vi^ that Xi ^ '='(^1 V^)- Suppose 
that rankPj^ = 1. Then we have either cardinahty of p equal to 1 or cardinahty of 
S equal to 1. Thus Ry^^ coincides with the Riesz measure px on ]0, +oo[ given by 

Px{^) = Y^l ^{u)u^-'du. (3.36) 

This implies that Xj = A and A > O which means that the result is true when 

rankPj^ = 1. Now, suppose that the claim holds for any i ^ p U S such that 
rankPj^ < k — 1, and let us show that it also holds for i such that rankP¿^ = 
k. Consider i E p U S such that rank7^¿^ = k. Then rankT'j^ = A; — 1 so that 
rankT'j^ < k — E Mi. Using the induction hypothesis, we have that Vj e M¿, 
= {A £ IR-) ^ £ Pi — 0; yi ^ Ij^} is in S(_7). Let R^. be a Riesz measure 
deñned as in (3.35) on the cone Vj^ for some í/j in and let C¿ = X! — {Pj ^ 

IR, j e /; Pj — O for j ^ h^}. Then from (3.14), the measure R^. = R(^., where 

is the convolution product, is concentrated on Vi^. Consider the sets 
Xi = {A, e IR, i e /; A, = O for j ^ 7¿^}, 
= {Pk £ IR, k & I; Pk ~ iT'ki for k e /¿^ and /S^ = O for A; ^ 



and 
where 



M{Xi) = {ak, k e /}, 



Xi for k = i 
^ for k e li^ 
O otherwise. 



Then it is easy to verify that M{Xi) e i3(i,'?/^i), where ipi G such that -01 (i) = 1, 
and -01 (j) = O Vj 7¿ i. Also we have Xj — M{Xi) = Xi ~ y ^ ^{h'^): where B{í,üj) 
is deñned by (3.26). 

Using the Laplace transforms, we obtain that 

Rxi = Rm{Xí) * R^._ñi- (3.37) 

Ai 2 
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Proceeding as in the proof of Theorem 3.2, and using (3.23), we get 

^ dimC' 

ZTT 2 r 

Rm{\í){v) = . / ^ (p{U.e^^)uf'~^duudCi, 

i- (Ai) -iTi+e^^ 



where U e 7^"*", U.e^^ = {ujk)j,kei ^-iid Ci is defined by (3.24). On the other hand, 
it is easy to verify that 



C/.e^i = ulEi + tí¿¿C¿ + II ll|¿ Ej- 

i-<3 



Henee, if we define 



(c„c,) ^ Q^(a,a)-Ell«.dl|i^.-, 

then we have 

Setting Un — y/v, we get 

dimC' 

2tT 2 /• + 00 p 

RmíxM = ^^^^ y'í^^í + v^C-, + Q\C,, Q))dCy^-'dv. 

This, using (3.36), becomes 

RMixd^) = ^'"^PkíI, Vi^^i + V^C, + Q'iQ, Ci))dQ),. (3.38) 
As for Xi — O, po is the Dirac measure at = O, we get 

Rmío){^)=R^^{^)=7i-^ í ^{Q\Ci,Ci))dCi. (3.39) 
Using (3.37) and (3.38), we obtain 

RM = 7r-"^pA,( / R^._n^ {^{vE, + + QXQ, Q) + Y))ydQ)),, (3.40) 

Denote by the set of C°° functions with compact support and consider the 
functions of the form 

ip{Z) = ipi{Zii)ip2{Zi^) 

where Z = {zij)ij^j G A, Zi^ e Ai^, (fii G C¿*(]R) and ip2 e C^{Ai^). Then by 
(3.39) and (3.40)', we have 

R^,{<p) = 7r-^pA,((^i) / Á._,.((^2(Q^(Q,a)+l^))yda 

2 

= PA,(V^l)(Áñi *i?.._2l)(¥'2) 
2 A.» 2 

= PA.(^i)i2xí(^2). (3.41) 

13 



For a suitable choicc of non-negative Lpi G C^(]R), we have Pa¿(v^i) > 0. If ip2 > O, 
then using (3.40) and the positivity of R^^, we get R^-_{ip2) = {p\-{(pi))~^ R^-{(p) > 0. 
Thus is positive and the induction hypothesis ensures that xi ^ '^{h 
Now, fix a non-negative (f2 such that R^..{ip2) is strictly positive. Then using again 
(3.40), we get PxX'^i) > O for any (pi > 0. Therefore px. is positive and we deduce 
that Aj > 0. If A¿ = O, then choosing aip ine^ such that — O, we get Xi £ S(i, ijj). 
To study the case Aj > O, we first observe that the map 

]0, +oo[xC' X {XeA,xu>0} 
is a diffeomorphism whose the inverse is given by 



1 



For a functions the functions (p e C^{A) of the form 



íPi{v)íp2{Ci)(ps{Y) {zn > 0), 
O {zu < 0), 



with (fii e C^{]0,+oo[), (p2 e C^{C') and (ps e C^{Ai^), by (3.40), we have that 

2 JC^ 



Since Aj > O, the positivity assumption of -R^. yields that R. is positive. This 

by the induction hypothesis imphes that Xi~ ^ ^'^{h '^)- Finally, choose a in 
such that 



uj{j) for Vj ^ i 
1 for j — i, 



where uj e e*. Then, as Xi ~ y ^ "^{h ^oi j G /.¿^, wc have that Xj = ^ if 

= o and Aj > if a;(j) = 1. As A¿ > O and n¿ = O, then A» > This means 
that for such a ■0, we have that Xí ^ '^{h '4')- Henee Xi ^ '^{'i) and Proposition 3.4 
is proved. □ 

Proof of Theorem 3.3 (<í=) Let % = ^ G S, where Xi ^ ^(0 let 

Rx= f[ Rx., (3-42) 
iepus 
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is the positive measure defined from Proposition 3.4. Then, for 9 &V* 

LrÁ^) = n Ln^Á&) 

iepuS 

- n A«(r^) 
= A,(ri)- 

(^) We have Ln^yd) = A^(^-i), then using the fact ^xA^~^) = Ax(^"^)> 

¿epu5 

Proposition 3.4, and putting Rx = H Rxi^ such that for 9 e P*, Lr^,{9) — 
Ag(ri), we get Xi e 5(0- Therefore ^ = ^ X¿ e S. □ 

FoUowing the terminology used in the papcr by Hassairi and Lajmi (2001) in the 
case of symmetric matrices, we cali the measures i?^, defined above in terms of their 
Laplace transforms, Riesz measures on the homogeneous cone. These measures are 
divided into two classes according to the position of x £ S. A class of measures 
which are absolutely continuous with respect to the Lebesgue measure on V and a 
class concentrated on the boundary &P of V. 

Proposition 3.5 Let x = {Aj, ¿ e /} G A". Then R^ is absolutely continuous if 
and only if Xi> i & I. In this case 

^^"^^^ ^ f^Ax+x(^)lp(^)^^, (3-43) 

where x = {—ni, i G /} and r-p(x) = tt^^— r(Aj — ^). 

iei 

Proof =>) We have V = ^ 7^"'".ei., then Xi defined by (3.30) is equal to x¿- 

* 

Prom (3.42), we have Rx ^ H Rxi- Writing Z = UU* e V, U e then 

Z = ^^Zi where Zi is in 7^"'".ei., and using the proof of Proposition 3.4, we have 
iei 

that for Zi = C/¿ = Mkjei 

RxÁdZi) = —^—-A^,{ZiMdZi), 

where A^'; = A^, and v{dZi) = vi,{dZi) = 1^%{Z) \{ du^j. 

^ ^ J ^ k 

m = 1 
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Then u{dZ) = A^{Z)dU, where 

1 + n • 

X = {-(^-^), j e 1} and dt/ = n n dukj. 

iepus i^j^k 

m = 1 

Using the fact that the mapping U G T^"*" ^ UU* EVisa diffcomorphism, we have 
that dU = 2-1-^1 A^(Z)(iZ, where x = {-^^, Vj e /}. As nj = l + + rij.), we 
have 

where x = {-nj, j G /} and Trix) = 2l^l ü r^+.e,(x.) = tt^ nr(A. - 

¿epu5 i€i 
Moreover, the condition A,¿ > i G / is easily deduced from Theorem 3.2. 

It suffices to verify that for % such that A¿ > i G /, the Laplace transform 

of the measure 

^;-^A^+^(Z)lp(Z)dZ 

is equal to A^(6'-^), Vó^ G V*. In fact, let 6» G "P*, thcn there exists T = (íj^) in 
7;+ such that e = T*T. Let Y = 7r(T)(Z), where tt is defined by (2.2). Then 
dZ ^detn-\T)dY and 

A^+^(z) = A^^^{n-\T)Y) = A^+^(r^y) = A^+^(r^)A^+^(y). 

Prom Andersson and Wojnar (2004), we have 

det7r-i(r) = ní.f'^'=A_^(0, 
iei 

where — x = {rii, i G /}. Writing Y — '^Y^, where y¡ is in 7^"^.ei. , then as 
rp(x) = 2l^l n r^+.e.,(x.) = 21^1 n L exp{-trr,}AÍ;(F,)z.i^(dr,) 

= |^exp-{try}A^+^(y)dy, 

we obtain that 

-1^/ exp{-tr(^Z)}A,+^(Z)dZ = A,(ri)-^ / exp -{trr}A,+^(r)c¿y 

rp(x) rp(x) 

= A,(ri). 

□ 
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4 Riesz exponential families 



In this section, we study the natural exponential family generated by a Riesz mea- 
sure. We first review some basic concepts concerning exponential families and their 
variance functions and introduce some notations. 
For a positive measure on A, we denote 

©(yu) = interior{^ e A*; L^{9) < 00} 

kf, = log 

where and /c^ are respectively the Laplace transform and the cumulant generating 
function of ¡i. 

The set M.{A) is now defined as the set of positive measures ¡j, such that ¡j, is not 
concentrated on an afñne hyperplane of A and ©(//) is not empty. For ¡i in A4{A), 
the set of probability 

F = = {P{e, n)dX = exp{-tr(^J5í) - k^}ii{dX)] 6 e ©(//)} 

To each ¡j, e Á4{E) and 9 G 6(/i), we associate the probability distribution on A 

P{e, ix){dX) = exp {{6, X) - k^{e)) fi{dX). 

The set 

F = F{f,) = {p{e,i2y, eeeii^)} 

is called the natural exponential family (NEF) generated by /x. We also say that 11 
is a basis of F. Note that a basis of F is by no means unique. If ¡i and u are in 
A4{A), then it is easy to check that F{fi) = F{u) if and only if there exist a & A 
and 6 G IR such that di'{X) — exp{{a,X) + h)d^{X). Therefore, if is in M{Á) 
and F = F{iJ,), then 

BF^{i^e M(A); F{v) = F} = {exp((a, X) + h)ii{dX)- (a, 6) G ^ x IR} 

is the set of basis of F. 

The function k^ is strictly convex and real analytic. Its first derivative k'^ defines a 

diffeomorphism between 0(/i) and its image Mp. Since k'^{9) = j XP{9, ¡i)(dX), 

Mp is called the domain of the means of F. The inverse function of k'^ is denoted 
by '?/^^ and setting P{m,F) = P{ipi^{m), /i) the probability of F with mean m, we 
have 

F = {P(m,F);m G Mp} , 
which is the parametrization of F by the mean. 

If ¡jl and u — exp((a, X) + b)/! are two basis of F, then for all 9 G D{u) — D{ij) — a, 

K{9) = k^{9 + a) + 6 (4.44) 

and for all m G Mp, 

— V'/ií^) ~ ^- (4.45) 
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Now the covariance operator of P(m, F) is denoted by VfÍitl) and the map defined 
from Mp into Ls{A) by m i — > Vp(m) = k'^{il)^{m)) is called the variancc function 
of the NEF F. It is easy proved that VF{m) = {ip'nim))^^ and an important feature 
of Vp is that it characterizes F in the foUowing sense: If F and F' are two NEFs 
such that Vir(m) and Vp/(m) coincide on a nonempty opcn sct of Mp fl Mp/, then 
F = F'. In particular, knowledge of the variance function gives knowledge of the 
NEF. 

Let íp{X) = S{X) + 7 be an afñnc transformation on A, where 6 G GL{A) and 
^ & A and let F be some NEF, generatcd by /x we denote by /Xj = ip * ¡j, the image 
measure of ¡j, by (p, then for all 6 e = 6* ^(©(^u)), 

k^^{e)^k,{6*{e)) + {e,^)- (4.46) 

The foUowing theorem gives a necessary and sufficient condition on x so that 
generates a natural exponential family. 

Theorem 4.1 Let x — {\, i & 1} be in S. Then the Riesz measure is in A4{A) 
if and only if Aj ^ O, for i & pU S . 

Proof (<í=) Suppose that x = {^i: i G /} is in S such that A¿ 7^ O V¿ e p U 5. 
We have Q{Rx) is not empty since it contains V*. We nccd to show that R^^ is not 
concentrated on an afñne hyperplane of A. Write X — Xi where Xi ^ ■^(0 (see 

* 

3.29). Then R^ = Y\. ^Xii sufñces to show that for any i e p U »S, R^^ is 

¿GpU5 

not concentrated on a afñne hyperplane of Ai^. In fact suppose that there exists 
i e p U íS such that R^^ is concentrated on a affine hyperplane H of Ai^. Then 
there exists G such that T^^ .e^ C H . On the other hand, there exist an element 
a G Ai^ and an hyperplane Hq of Ai^ such that H = a + Hq. Write 

Ai^ — Hq + IRX, 

where X ^ Hq. Let (efj)i<fe<„;^. be a basis of Aij, G /¿^ x 7j^. 
As Ai^ — Y\. '^iji write 

^= E E 

l,jeli^ l<k<nij 

where for /, j G /¿^ and 1 < /c < n^j, is a real number. Thus, we can write 

A^^Ho+Y. E ^4- 

Since the dimensión of Hq is equal to dim(^j^) — 1, then there exist l,j G /j^, and 
1 < A; < riij such that 

Ai^ = Ho + JReí. 
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Let US consider the vectors 

It is clear that Ai and A2 are in .e^ C ií. Using the fact that for ¿ e p U <S, 
Aj 7^ O, we have necessarily 7^ O, and we get A2 — Ai = efj- which is an element 
of Hq. This is in contradiction with the fact that X ^ Hq. Thus for any i ^ pU S, 
R^. is not concentrated on a affine hyperplane of Ai^ and Thereoem 4.1 is proved. 

Suppose that Xi = 0,\/i E pU S. As the support of is P, thus R^ is not an 
element of Ai (^). □ 

Next, we give the variance function of the Riesz exponential family F{R^) generated 
by R-^. For X and K in A, we define the quadratic representation P{X) by 

P{X)K = X{KX)- 

It is symmetric, since we have {P{X)K, L) — {P{X)L, K). 

Theorem 4.2 For any m eV, 

For the proof of this theorem we were led to estabhsh the foUowing intermediary 
result 

Lemma 4.3 The map 

(p:X^ X-^ -.V ^ B 

Y ^ tr{X-^Y) 

is differentiable and its differential is 

ip'{X){K) = {X-^YiK) = -P{X-'){K), 

that is 

{X-')\K){Y) = -tr{{X-\KX-'))Y). 

Proof We have 

liml((X + íX)-i-X-i)(y) = \\mUr{{X-rtK)-\{X-{X-rtK))X-^)Y) 

= -\mv^ii{{X + tK)-^{KX-^)Y) 
= -{X-\KX-')){Y) 

We now show that 

{X + K)-^ - X-^ + X-^{KX-^) = o{K) 
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\\{X + K)-^ -X-^ + X-^{KX-^)\\ = \\ {{X + Ky^{X - {X + K) + {X + K){X-^K))X 

= II ((X + K)-\X -iX + K) + {X + K){X-^K))X 
= \\{{X + K)-\K{X-'K))X-^)\\ 

< II (x + ií)-^ lili lili íí^ II . 

Then 

{X-\K)y = -X-\KX-^) = -P{X-^){K). 

□ 

Proof of Theorem 4.2 For 9 e V*, we have Lr^{9) = A^(^-^). Then using (3.10), 

we get 



Therefore 



iei 

and from Lemma 4.3, we get 

It is easy to see that (O^i)"^ — {0^^)^i and {0^i)~^ — {9~^)^i, then 

For m eV, such that m — k'ji^{9), we have 

iei 



then 



-m. 



¿^7 \" /t^ . ■ 

i '^i 

Therefore 



□ 



Note that, in the particular case of the Wishart NEF, that is when x — {K^i € I}, 
we have 
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